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We present a Bayesian estimation analysis for a particular trace gas detection technique with 
species separation provided by differential diffusion. The proposed method collects a sample con- 
taining multiple gas species into a common volume, and then allows it to diffuse across a linear 
array of optical absorption detectors, using, for example, high-finesse Fabry-Perot cavities. The 
estimation procedure assumes that all gas parameters (e.g. diffusion constants, optical cross sec- 
tions) are known except for the number population of each species, which are determined from the 
time-of-flight absorption profiles in each detector. 



Trace gas detection techniques have wide-ranging ap- 
plications including explosives detection [T] and medical 
diagnostics [2^ . Individual sensors within a given system 
should have adequate sensitivities for measuring physical 
properties of the gas, but must also demonstrate the se- 
lectivity necessary to distinguish between the component 
gas species. For example, if the optical absorption by a 
uniformly mixed sample of two gas species is measured 
at a single optical frequency, one can only determine the 
(scaled) sum of the two populations, but not the indi- 
vidual populations. There are several approaches for ob- 
taining selectivity, including the use of multiple probe 
frequencies along with a knowledge of the different gas 
spectra. This can be done with sequential measurements 
with a tunable laser ^ or simultaneously, as in recent 
work with femtosecond lasers [2]. Alternatively, one can 
spatially separate the species as in, for example, the in- 
dustrial workhorse methods of mass spectrometry and 
gas chromatography. 

Given a method that is both sensitive and selective, 
the task remains to estimate the populations of each gas 
within the original sample from the raw measurement 
data. Bayesian estimation provides the most natural 
framework to update one's state of knowledge accord- 
ing to the acquired data. In this approach, we start with 
a prior number distribution for each species, then iter- 
atively update the entire joint probability distribution, 
using the measurement results and Bayes' rule [3J. The 
widths of this joint distribution will narrow in time due 
to this procedure, reflecting the acquisition of knowledge. 
As long as the assumptions going into the experimental 
model are accurate, the Bayesian approach will be opti- 
mal, in contrast to more common heuristic approaches. 
Bayesian methodology is becoming increasingly relevant 
in many areas of physics and technology [U [S] as ex- 
periments reach such levels of sophistication that crude 
estimation techniques severely limit performance. For ex- 
ample, in single-particle tracking experiments, Bayesian 
estimation can be used to extract the particle type based 
on the estimation of its diffusion constant [4 . 

In this paper, we propose a physical device capable 
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of selectively measuring individual trace gas concentra- 
tions from a mixture containing multiple gas species, 
based on characteristic diffusion times. We also present 
the Bayesian estimation formalism for mapping the de- 
vice's raw measurement output to the individual trace 
gas populations. The essential configuration of the dif- 
fusion spectrometer we propose is displayed in fig. lA. 
First, a sample containing a number of gas species is 
collected within the initial common volume and then re- 
leased to diffuse into the sensor chamber, which is filled 
with a buffer gas and lined with a linear array of detec- 
tors. We consider explicitly a series of high-finesse optical 
cavities which measure optical absorption at a single fre- 
quency. However, the Bayesian method presented here 
can be easily extended to different trace gas detectors 
beyond the optical cavity, with varying forms of space- 
or frequency-based selectivity. The difference in diffu- 
sion constants causes the different species of atoms to 
transiently fill the chamber at differing rates, allowing 
for a selectivity not possible in a configuration where all 
of the species are mixed homogcnously. Fig. IB dis- 
plays the measurement output for each cavity along with 
the contributions from each gas at each point in time. 
Note that even a device with only one cavity displays se- 
lectivity due to differential diffusion, however increasing 
the cavity number is beneficial both for the extra infor- 
mation, and for allowing common-mode cancellation of 
technical cavity noise (which is a common motivation for 
gradiometer configurations) . 

More specifically, an initial volume of gas containing 
Ni particles of gas i is released into the chamber. We 
assume each gas has a distinct diffusion constant Di (in 
units of m^/sec) and optical scattering cross section di 
(in units of m^), and that these quantities are known for 
each potential gas in our library. We assume that the 
gases uniformly fill the transverse direction, such that 
the distribution is only given along the direction of prop- 
agation as 

The number density at x is given by pi{x,t) ~ 
NiPi{x,t)/A where A is the cross-sectional area of the 
tube. The optical absorption coefficient is defined as 
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FIG. 1: (A) A gas mixture is released from a loading chamber 
and diffuses across two cavity modes, each of which gener- 
ates an absorption signal. (B) The noisy plots show the raw 
absorption output from the two cavity measurements. The 
dashed (red) curves denote the contribution from gas-1 and 
the dashed-dotted (blue) curves denote the contribution from 
gas-2. The set of curves which peak first are for cavity-1 and 
the set which peak second are for cavity-2. The chamber has a 
square uniformly filled transverse cross section 1 cm per side. 
Each cavity mode has a length 8 mm, and a transverse width 
(along x) of 0.6 mm. The cavity modes are located along x 
at a distance of 1 and 2 mm from the origin. The gas is col- 
lected into a collection volume of size 10 mm^ (0.1 mm width 
along x), where the initial concentrations are set at 1 and 2 
ppb (parts per billion). The optical absorption cross-section 
is the same for each gas and set to at = 10~^^ cm^ which 
is a typical value for the absorption of infrared radiation by 
molecules such as ethane {C2H6). The diffusion constant for 
gas-1 is Di = 0.044 cm^/s (the approximate value for ethane 
at room temperature and atmospheric pressure) and gas-2 is 
set a factor of two larger. The white measurement noise cor- 
responds to a (jm ~ 10^^^ cin~^ /VHz absorption sensitivity. 

ai{x,t) — pi{x,t)ai (in units of 1/m). The width of the 
distribution at time t is given by erf (t) — erf (0) -f 2Dit. 
We place a number of sensors in a hnear array along x, 
and the measurement output for a sensor j at location 
Xj is given by 

Vjit) = + (2) 

i 

where the absorption at Xj due to gas i is ai{xj,t) = 
Cj^i{t)Ni, Cj^i{t) = pi{xj,t)ai/A, and £,j{t) is a white 
noise term (of variance 1/dt). We assume that each 
detector has the same absorption sensitivity ctm, which 
for high-finesse optical Fabry- Perot cavities can approach 
10-14 cm-i/VSz [6J. 

In full generality, the Bayesian approach would use 
each successive measurement result, the assumed model, 
and B ayes' rule to update the collective distribution that 
describes our state of knowledge about the gas number 



populations. This description can be simplified if we as- 
sume that the distribution begins and remains Gaussian. 
In this case, we do not have to evolve the entire distribu- 
tion, but merely a small set of moments, or more specifi- 
cally, the number estimates TVg (the distribution means), 
and the covariance matrix, S(t) (the distribution widths) 
[3] . The diagonal elements of the covariance matrix self- 
consistently represent our uncertainty in the populations 
and the average error of our estimates, {{Ne,i — ^i)^), 
which the Bayesian estimator serves to minimize. 

Technically, this Gaussian approach is not always ideal. 
First, the Gaussian tail allows for the possibility of a non- 
physical negative population. Second, for many trace gas 
scenarios the distribution is naturally bimodal with one 
peak at zero (no gas present) and another peak at a non- 
zero value (at the typical concentration level when the gas 
is present). Nevertheless, the Gaussian approximation is 
illustrative and simplifies the numerical treatment. To 
indicate that a certain level of gas is 'detectable' in this 
simplified framework, we require that the width of the 
distribution be less than the typical concentration value 
when the gas is present. 

The covariance matrix is initialized at S(0) = Sq which 
is a diagonal matrix with elements corresponding to the 
prior level of number uncertainty for each species. Ac- 
cording to the Bayesian procedure with the Gaussian as- 
sumption, the equation of motion for the covariance ma- 
trix, or Riccati equation, is the deterministic nonlinear 
matrix differential equation 

dS(t) _ S(<)C^(t)C(t)S(t) 



This solution of this equation can be written [7] 

When the initial spatial gas distributions are approxi- 
mated as delta functions (erf (0) = 0), the integral can be 
represented in terms of the exponential-integral function. 

This treatment results in the following estimator which 
maps the noisy measurement results to the best-guess 
estimates for the population of each species in the sample 
volume: 

= Ko{t) [m - cm^it)) . (5) 

Here y{t) — C{t)Ne{t) represents the 'innovation' or the 
extent to which the measurement results deviate from our 
expectation. The matrix Ko{t) = Y.{t)C'^ {t) / a\j repre- 
sents the 'observer gain' which weights the importance of 
the new measurement results [3]. Note that the covari- 
ance matrix not only reflects the expected performance 
but also tells us how to tailor this observer gain in time. 
Namely, the more information that is acquired, the less 
relative importance we give to the latest data point. 

The Bayesian estimator performance is shown in fig. 
2. Here we introduce two species into the chamber 
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FIG. 2: (A) The Riccati solution of eq. ^ for the scenario and parameters of fig. 1, but with cta/ — 10^^^ cm~^/\/Hz. The 
number uncertainty for gas-1 is plotted as a dashed (red) line, and for gas-2 as a dashed-dotted (blue) line. The horizontal 
line denotes a target 100 ppt (parts per trillion) detection level. If there were no uncertainty in the population of gas-2, the 
estimation performance for gas-1 would improve to that shown by the dotted (red) line. If the diffusion constants were equal 
or the gases were immobile, the estimation performance would hardly improve upon the initial uncertainty. (B) For a single 
trial, the population estimates as a function of time are shown. The estimate for each is surrounded by curves separated by 
the Riccati solution uncertainty. The true value for each gas is given by the horizontal lines (corresponding to 1 and 2 ppb) , to 
which the estimates converge at long times. 



with equal optical scattering cross-sections (cti = (12), 
but with diffusion constants differing by a factor of two 
(Z?2 = 2Z?i) and initially unknown concentrations. Note 
that we have chosen a conservative value for ctm = 
10^^^ cm^^/VHz, and thus significantly improved per- 
formance should be possible before reaching a typical 
photon shot noise limit. We assume prior uncertainty in 
the number of each gas at 10 ppb, and as time increases, 
we determine the number in each species with increasing 
confidence. The uncertainty variances eventually level off 
as the gas dissipates throughout the chamber and reaches 
a flat steady-state distribution. We recognize that in the 
absence of uncertainty in gas-2, the measurement of gas- 
1 improves as denoted by the dotted line. However, as 
the uncertainty of gas-2 is increased, reflecting the pos- 
sibility of a gas sample with trace gas-1 and large back- 
ground concentration of gas-2, the performance of the 
gas-1 measurement remains near the level indicated by 
the dashed line. This indicates that the method will work 



well to identify trace components despite the possibility 
of a large background concentration. 

In conclusion, we have used Bayesian estimation to 
show that an array of absorption detectors can be a 
highly selective trace gas detector without the need for 
frequency tunability. While we have presented a simpli- 
fied Gaussian treatment for numerical purposes, a full 
Bayesian treatment with bimodal distributions and 3-D 
diffusion requires a simple extension of these concepts. 
For this diffusion selectivity approach, increased pre- 
cision in the library of known quantities for each gas 
species (e.g., Di, ai) under various environmental con- 
ditions would improve performance. With the technical 
challenges of existing selectivity methods in trace gas de- 
tection including expense, portability, and bandwidth, 
our diffusion model with Bayesian estimation presents a 
practical alternative. 
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